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Abstract 

I We show simple methods how to separate pure CP violating effect from matter 

J> ' effect in long baseline neutrino oscillation experiments with three generations of 

neutrinos. We give compact formulae for neutrino oscillation probabilities assuming 
one of the three neutrino masses (presumably u T mass) to be much larger than the 
£f~) | other masses and the effective mass due to matter effect. Two methods are shown: 

One is to observe envelopes of the curves of oscillation probabilities as functions of 
. neutrino energy; a merit of this method is that only a single detector is enough to 

determine the presence of CP violation. The other is to compare experiments with 
O 1 at least two different baseline lengths; this has a merit that it needs only narrow 

q^. energy range of oscillation data. 
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1 Introduction 

X 

The CP violation has been observed only in the hadron sector, and it is very hard for us 
to understand where the CP violation originates from. If we observe CP violation in the 
lepton sector through the neutrino oscillation experiments, we will be given an invaluable 
key to study the origin of CP violation and to go beyond the Standard Model. 

The neutrino oscillation search is a powerful experiment which can examine masses 
and/or mixing angles of the neutrinos. The several underground experiments, in fact, have 
shown lack of the solar neutrinos pi |3j] and anomaly in the atmospheric neutrinos 
|| H, strongly indicating the neutrino oscillation |10|, [n], [L2| . The solar neutrino deficit 
implies a mass difference of 10~ 5 ~ 10 _4 eV 2 , while the atmospheric neutrino anomaly 
suggests a mass difference around 10~ 3 ~ 10~ 2 eV 2 [[TO], [□], [PJ]. 

The latter encourages us to make long baseline neutrino oscillation experiments. Re- 
cently such experiments are planned and will be operated in the near future [O, [TJJ. It 
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1 Some experiments have not observed the atmospheric neutrino anomaly |j, 
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is now desirable to examine whether there is a chance to observe not only the neutrino 
oscillation but also the CP or T violation by long baseline experiments [[15|, |TjJ. Two of 
the present authors have studied how large T violation we may be seen in long baseline 
experiments [|16|], but they have not answered the question how the CP violation is distin- 
guished from the matter effect, in case both the solar neutrino deficit and the atmospheric 
neutrino anomaly are attributed to the neutrino oscillation (In this case the matter effect 
is expected to give a fake CP- violation effect comparable to pure CP- violation effect). 

In this paper we will answer this question. In sec. 2 we briefly review neutrino oscilla- 
tion for 3 generations, and give very compact formulae describing neutrino oscillation in 
the presence of matter (The detailed derivation of the formulae is given in the Appendix). 
In sec. 3 we show two methods to distinguish pure CP violation from matter effect. In 
sec. 4 we summarize our work and give discussions. 



Compact Formulae for the Neutrino Oscillation Prob- 
abilities 



2.1 Brief Review and Parameterization 



Let us briefly review CP violation in neutrino oscillation [17|, [L8], [iij to clarify our notation. 

We assume three generations of neutrinos which have mass eigenvalues m,i(i = 1,2, 3) 
and mixing matrix relating the flavor eigenstates v a [a = e, /x, r) and the mass eigen- 
states in the vacuum v'^i = 1, 2, 3) as 



(1) 



We parametrize [EC, 21, 22] with the Gell-Mann matrices Aj's as 



Jj(0) _ e iiM 7 -p e i</>A 5e ia;A2 

1 \ / 1 

Si> I I 1 

-s^p Cip 



e" 




AS 



(2) 



where = cosip, = sin0, etc. 

The evolution equation for the flavor eigenstate vector in the vacuum is 



.du 

dx 



-U^di ag ( Pl , P2 ,p 3 )U^u 

-Pi + ^ (0) diag(0, 5m 2 21 , 5m 2 31 )U^} 



v. 



(3) 
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where p^s are the momenta, E is the energy and Smfj = mf — m 2 . Neglecting the term 
Pi which gives an irrelevant overall phase, we have 



i^ = ^Wdiag(0 l tf<,*O^°)ti/. 



Similarly the evolution equation in matter is expressed as 



i— = Hu, 

ax 



where 



1 



(4) 

(5) 
(6) 



with a unitary mixing matrix U and the effective mass squared fif's (i = 1,2,3). The 
matrix U and the masses /ij's are determined by 



(7) 
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Here 



a = 2V2G F n e E 
= 7.56 x 10" 5 eV 2 



P 



E 



(8) 



g cm 3 GeV 

where n e is the electron density and p is the matter density. The solution of eq.(|5]) is then 

v{x) = S(x)u(0) (9) 

with 



S = Te 



-if X dsH(s) 



(10) 



(T being the symbol for time ordering), giving the oscillation probability for v a — > 

vp(a, P = e, fx, t) at distance L as 

P{v a -^v (i -L) = \Sfs a (L)\ 2 . (11) 

The oscillation probability for the antineutrinos P{p a — > %) is obtained by replacing 
a — >■ -a and U — > U*(i.e. 5 — > -5) in eq.(|TTD. 

We assume in the following the matter density is independent of space and time for 
simplicity, and have 



S(x) 



-iHx 



(12) 
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2.2 Approximation of the Oscillation Probability 

If we attribute both the solar neutrino deficit and the atmospheric neutrino anomaly to 
the neutrino oscillation with MSW solution for the solar neutrinos, we find most plausible 
solutions to satisfy 8rrv^ x <C 8rrv^ l and a <C 8rrv^ 1 |H| [T1J. In the following we assume 
a, 5m\i <C 3 m li Q- This case is also interesting when we consider the long baseline 



neutrino oscillation experiments to be done in the near future |13| , p~4 
Decomposing H = H + Hi with 



/ 



IE 



tf(0)t 



and 



2E 



( 



V 



5m 2 2l 



5m 2 31 



(13) 



(14) 











we treat Hi as a perturbation and calculate eq. (|T2|) up to the first order in a and 5m| x . 
Defining Q(x) and Hi(x) as 



and 

we have 
and 

which give the solution^ 



Hi(x) = 



jH x 



S(x) 



jH x Hie -iH x 



i — = Hi(x)n(x) 

ax 



fl{x) 



fi(0) = 1, 

: 1 - i pdsH^s). 
Jo 



(15) 
(16) 

(17) 

(18) 



(19) 



We note the approximation (|T9| ) requires 



ax 
2E 



<C 1 and 



5m\iX 



2E 



< 1. 



The equations (|T5|) and ([TO]) give 

S(x) ~ e-' lHoX 



+ e 



" iHox (-i) r AsHAs). 
Jo 



(20) 



(21) 



2 For the case Sm^i <C a -C om^ see ref.|16||. 

3 We note the eq.(|l^) is correct for a case where the matter density depends on 
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We then obtain the oscillation probabilities P(f M — > v c ), P(^ —> v^) and P(f M — > u T ) in 
the lowest order approximation as 



P(i/ jU -> z/ c ; L) = 4 sin 
(5ml L . 5m%i L 



5m%i L 



2 " ""31^ 2 2 2 I i 

4 £ c ^ s </, i 1 



5m 2 31 



2(1 - 2sJ) 



+ 2 " „ 3 i sin „ 3 i c|g^ 



2 5rr4 



<*31 



- 4 



2P 2E 

SrriLL 2 Jm^L 2 
2E sm 4 ^ s^c^c^c^ 



— — s^(l - 2s ) + — f L S^(-S0Sv S ^ + C <5 C V C ^ 

0777.Q1 712 



31 



(22) 



PW - ^) = 1 + 4 sm 2 — ^—44 { (44 - 1) + — . 2 S j(i - 244) 



+ 2 



4P 

(Jm^ L . L 5ml 



2E 



Smix 



4(244 - 1) 



+ 



^ m 21 i 2 2 .2 , 2 . 2 



and 



• " ) 4 sin 2 ^l^cic 2 , 2 | I • I. 



4P 



+ 2 ^ sm JJ: c^s^ 



IE 
Sm 



IE 



2 ^4 c V'4 s V' 



31 

+ 4 2 ^ sm ssc^c^s^s^. 



(23) 



(24) 



(Detailed derivation is presented in the Appendix). Recalling that P{y a — > i/g) is obtained 
from P{y a — > i/g) by the replacements a — > —a and 5 — > —5, we have 



P(z/ M -> u c ; L) - P(z/ M -> u e ; L) 

APi(^ -> I/ e ) + AP 2 (Z/ M -> I/ e ) + AP 3 (Z/ M -> I/ e ) 



with 



and 



AP 1 (u^u c ) = 
AP 3 (z/ M ^z/ c ) = 



2 HP VvM 1 zs <t» 



16 =- sm 

5m§! 4P 



-4|§sin ^444(1-24), 



(JmoiL 2 5ml-,L 2 
' 2 ^ sm ^ ssc^c^s^s^. 



(25) 

(26) 
(27) 

(28) 
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Similarly we obtain 



a 

16 



5m 2 31 



sm 



2 5 m 2 31 L 
AE 



1 5m% L 5m% L 
"- sin 



A 2E 



2E 



2 2 2 

C <f> S (/> S ip 



2r 2 s 2 



(29) 



and 



+ 



AP(u, 
-32 



a 



<5m 2 



z/ T . 



sm 



31 



2 5m|iL 
4P 



1 5ml L (5mq, L 

— sm — — 

A 2E 2E 



5m%iL 9 5m 2 .-, L 
' 21 sin 2 31 



2E 



AE 



4 2 2 2 

C (f> s 4> c ip s ip 



(30) 



Here we make some comments. 



1. P[y a — > vpYs and AP(u a — > i^)'s depend on L and E as functions of L/P apart 
from the matter effect factor a (= 2y2Gpn e E). 

2. At least four experimental data are necessary to determine the function AP(z/ M — > 
z/ e ), since it has four unknown factors: 5m 2 1; 5m 2 1? c^s^s^l— 2s 2 ) and Sgc^s^c^s^c^s^ 
In order to determine all the mixing angles and the CP violating phase, we need to 



observe P{v„ 



z/ M ) and P(z/ M 



Up) in addition. 



3. AP(i/ p 



u^) is independent of 5 and consists only of matter effect term. 



3 Separation of Pure CP Violating Effect from the 
Matter Effect 

Next we investigate how we can divide AP(z/ M — ► u e ) into a pure CP-violation part and a 
matter effect part The terms AP 1 (z/ /i — > z/ e ) and AF 2 (^ — > fe), which are proportional 
to "a", are due to effect of the matter along the path. The term AF 3 (i/ (J — > u e ), which is 
proportional to s§, is due to the pure CP violation (We simply call AP^v^ — > u c ) as APj 
hereafter). In the following we introduce two methods to separate the pure CP violating 
effect AP3 from the matter effect APi + AP2. 



3.1 Observation of Envelope Patterns 

One method is to observe the pattern of the envelope of AP, and to separate AP3 from 
it. Considering the energy dependence of a(oc E), we see that APi/L, AP 2 /L and AP 3 
depend on a variable L/E alone. The dependences of them on the variable L/E, however, 

4 It is straightforward to extend the following arguments to other processes like — > v T . We present 
the cases of — > v e and — > D e as examples. 



6 



are different from each other as seen in Fig. Each of them oscillates with common zeros 
at L/E = 2im/5m\ 1 (n = 0, 1, 2, ■ • •) and has its characteristic envelope. The envelope of 
APi/L decreases monotonously. That of AP 2 /L is flat. That of AP 3 increases linearly. It 
is thus possible to separate these three functions and determine CP violating effect AP 3 
by measuring the probability AP over wide energy range in the long baseline neutrino 
oscillation experiments. This method has a merit that we can determine the pure CP 
violating effect with a single detector. 

In Fig. ^| we give the probabilities Piy^ — > u e ) and P{v^ — > v e ) for a set of typical 
parameters which are consistent with the solar and atmospheric neutrino experiments 
|TTJ: 8m 2 21 = 10- 4 eV 2 ,5m^ = HT 2 eV 2 ,s^ = 1^/2, s = V^l and = 1/2. We see 



the effect of pure CP violation in Fig. 0(a), since we find that the curve AP has the 
envelope characteristic of AP 3 . We show in Fig. ||the same probabilities as Fig. ||(a) but 
as functions of E to see the energy dependence more directly. 

We comment that the envelope behavior of AP can be understood rather simply as 



follows: The term AP3 is proportional to [23, 16 



/ = sin A 2 i + sin A 32 + sin A 13 

. A 2 i . A 32 . A13 
= — 4 sin sin sin , (31) 



where 



A = 32 
3 2E ( J 

Since we are interested in the first several peaks of /, we have 

A 31 ,A 32 ~0(1). (33) 

On the other hand we have 

A21 < 1, (34) 
because A 2 i A31, A 32 . Taking into account 

A 2 i + A 32 + A 13 = (35) 

and eqs. (|33|) and fl34|), we obtain 

AP 3 cx/~2A 21 sin 2 ^. (36) 

This shows AP 3 has a linearly increasing envelope A 21 oc L/E. On the other hand, the 
envelopes of APi and AP 2 do not increase with L/E for fixed L, and it makes AP 3 
dominant in AP for large L/E. 
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(a) Matter effect term APi(z/ M -> v c ) divided by L for cjsjsj(l _ 2 4) > °- The 
envelope decreases monotonously with L/E. 





(c) CP-violation effect term AP 3 (z/ At — > z/ c ) for ssc^s^c^s^c^s^ > 0. The envelope 
increases linearly with L/E. 

Figure 1: The oscillation behaviors of the APi, AP 2 and AP 3 . 
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(a) The oscillation probabilities as functions of L/E for 5 = tt/2. 
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(b) The oscillation probabilities as functions of L/E for 5 = 0. 



Figure 2: The oscillation probabilities for S = rt/2 (Fig. 0(a)) and 5 = (Fig. 0(b)). 
P{ v n ~^ Z7 e),-P(^ —>■ ^e) and AP(^ — > u e ) are given by a broken line, a dotted line 
and a solid line, respectively. Here p = 3gcm~ 3 and L = 250 km (the distance between 
KEK and Super-Kamiokande) are taken. Other parameters are fixed at the following 



values which are consistent with the solar and atmospheric neutrino experiments UTT 



Sm^ = 10~ 4 eV 2 ,Sm 2 31 = lQ- 2 eV 2 ,s^ = l/y/2, = ^fOA and Suj = 1/2. 
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Energy / GeV 



Figure 3: The oscillation probabilities as functions of E. Parameters are taken the same 
as in Fig. 0(a). 
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3.2 Comparison of Experiments with Different L's 



The other method is to separate the pure CP violating effect by comparison of experiments 
with two different L's. Suppose that two experiments, one with L = L\ and the other 
L = L2, are available. We observe two probabilities Piy^ — > u e ; Li) at energy E\ and 
P{ u ii ~ * v e ; L2) at energy E2 with L\/E\ = L2/E2. Recalling that P[y^ — > v e ',L) is a 
function of L/E apart from the matter effect factor a(oc E), we see that the difference 

{Pfa -> i/ e ; Li) - P(z^ -> i/ e ; ^2)} Ll/jBl=i2/£2 (37) 



is due only to terms proportional to 
(APi + AP 2 ) from AP(^ -> i/ e ) asQ 



"a". 



We obtain AP 3 by subtracting these terms 



AP 3 K 
AP(V„ 



24; Li) 
• Li) - 



2Li 



AP(z/ M -> i/ e ;L x ) 



L2 — L\ 
L 2 — Li 



P(v» -> Li)} 



L/E=const. 



(3? 



{AP(z/ M -> z/ c ; La) - AP(^ -> z/ e ; L^} 



(39) 



L/_B=const. 



This method has a merit that it does not need to observe the envelope nor many oscillation 
bumps in the low energy range. 

In Fig. [| we compare Piy^ — > v e ) for L = 250km (KEK/Super-Kamiokande exper- 
iment) with that for L = 730km (Minos experiment) in a case with the same neutrino 
masses and mixing angles as those in Fig. 0(a) (or Fig. [3]) . We see their difference, 
consisting only of the matter effect, has the same shape as the solid line in Fig. 0(b) up 
to a overall constant. We also show the pure CP violating effect obtained by the two 
probabilities with eq.(B3). This curve has a linearly increasing envelope as seen in Fig. 



c . 



4 Summary and Discussions 

We have given very simple formulae for the transition probabilities of neutrinos in long 
baseline experiments. They have taken into account not only the CP-violation effect 
but also the matter effect, and are applicable to such interesting parameter regions that 
can explain both the atmospheric neutrino anomaly and the solar neutrino deficit by the 
neutrino oscillation. 

We have shown with the aid of these formulae two methods to distinguish pure CP 
violation from matter effect. The dependence of pure CP-violation effect on the energy 
E and the distance L is different from that of matter effect: The former depends on L/E 
alone and has a form f(L/E), while the latter has a form L x g(L/E) = E x g(L/E). 
One method to distinguish is to observe closely the energy dependence of the difference 
P{ u n ~* v c\ L) — P{p^ — > u e ; L) including the envelope of oscillation bumps. The other is 

5 Note that the eq.(^8|) does not require P{D I1 — * v c ; L 2 ). 
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Figure 4: The oscillation probabilities P(V M — > f e )' s f° r KEK/Super-Kamiokande exper- 
iment with L = 250km (broken line) and those for Minos experiment with L = 730km 
(dotted line). Masses and mixing angles are the same as in Fig. ||j(a). Their difference, 
which consists only of matter effect, is shown by a dot-dashed line. The pure CP violating 



effect in KEK/Super-Kamiokande experiment determined by eq. (|38|) is drawn by a solid 
line. 
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to compare results from two different distances L\ and L 2 with L\jE\ = L2/E2 and then 
to subtract the matter effect by eq.(p8|) or eq.([39|). 



Each method has both its merits and demerits. The first one has a merit that we need 
experiments with only a single detector. A merit of the second is that we do not need 
wide range of energy (many bumps) to survey the neutrino oscillation. 

It is desirable to make long baseline neutrino oscillation experiments with high inten- 
sity neutrino flux, and to study CP violation in the lepton sector experimentally. 

A Derivation of the Oscillation Probabilities 

Here we present the derivation of eq.(^) ~ eq.(^4|) with use of eq.(^), and show how 
well this approximation works. Let us set S(x) = Sq(x) + Si(x), defining 

S (x) = e~ iHoX : (40) 
S 1 (x) = e~ iHoX {-i) ["dsHtis). (41) 

We see 







= ^ + ^*(e-*-l) (42) 



and 



S 1 (x) fia = ~i£ds [e-^^^e"^ 

= -iU^U^m^fU^Tix)^ (43) 



where 



Si 3 S j3 ■ xe~'^E- 



,2 \-i 



Using 



+ (l-S i3 )(l-S j3 )-x. (44) 



U^m-M? = 7^ {diag(0, 5m 2 21 , 0) + U^dmg(a, 0, 

= ^M* + ^M> (45) 
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Figure 5: Exact and approximated values of P{p ia 
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for L = 250 km (Fig. |(a)) and 



those for L = 730 km (Fig. 0(b)). Exact values and approximated ones are shown by 
solid lines and white circles, respectively. The parameters s^, s^, s u , 5 and p are taken the 
same as in Fig. @(a). 
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and 

we obtain 
with 

iT(x)p a -- 



k=i 



S{x)/3a = 5p a + lT{x)p a 



(46) 
(47) 
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13 
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13 
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We then obtain the oscillation probability in the lowest order approximation as 



P(v a -> z^; L) = ^(L)^ 
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a 7 . bm^ x L 



-2 — sm 
277 2E 



U, 



(o) 

a 3 



17 



(o) 



13 



+ 4 



U, 



(0) 



+ 2- 



33 

~2E 

a 



U, 



(o) 



a 3 



Sm^L 



sm 



sm 



8m% 



- 4- 



'31 

, 5wioi 7 



277 

5 a i5/3i 



477; 
gmfi 



1-4 



(0) 
13 



+ 2 



(^al + 8(3l) 
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(0) 
13 



- 5 a i - 8 I 



pi 



277 



' oi*^Im (U%*U$U3U%*) 



477 



(49) 



Substituting eq.(^j) in eq. (^9|) we finally obtain eq.(|22|) ~ eq. 

Figure [5] shows how well this approximation works for KEK/Super-Kamiokande ex- 
periment and also for Minos experiment with the same masses, mixing angles and CP 
violating phase as in Fig. 0(a). Our approximation requires (see eq. (^0|)) 



— = °- 420 ( ~\ 1 ( P 3 

277 V 730 km/ I 3 gem" 3 



< 1 



and 



Sm^L 



0.185 



(5m| 1 /10- 4 eV 2 )(7/730km) 



« 1, 



(50) 



(51) 



277 77/GeV 

which is marginally satisfied for L = 730km. We see that even in this case eq.(|4"9"D gives 
good approximation. 
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